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Abstract. The rapidity-even dipolar flow vi associated with dipole asymmetry in the initial 
geometry is measured over a broad range in transverse momentum 0.5 GeV< pt <9 GeV, 
and centrality (0-50)% in Pb-Pb collisions at ^/sT^ = 2.76 TeV, recorded by the ATLAS 
experiment at the LHC. The vi coefficient is determined via a two-component fit of the first order 
Fourier coefficient v\,\ = (cos A0) of two-particle correlations in azimuthal angle A(f> — (f>a — <f>b 
as a function of p T and p T . This fit is motivated by the finding that the pt dependence 
of ui,i(pxiPt) data are consistent with the combined contributions from a rapidity-even vi 
and global momentum conservation. The magnitude of the extracted momentum conservation 
component suggests that the system conserving momentum involves only a subset of the event 
(spanning about 3 units in rj in central collisions). The extracted vi is observed to cross zero 
at pt ~ 1.0 GeV, reaches a maximum at 4-5 GeV with a value comparable to that for i>3, and 
decreases at higher pt- Interestingly, the magnitude of Vi at high pt exceeds the value of the V3 
in all centrality interval and exceeds the value of V2 in central collisions. This behavior suggests 
that the path-length dependence of energy loss and initial dipole asymmetry from fluctuations 
corroborate to produce a large dipolar anisotropy for high pt hadrons, making the vi a valuable 
probe for studying the jet quenching phenomena. 



Recently, the measurements of harmonic flow coefficients v n have provided valuable insights 
into the properties of the hot and dense matter created in heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). These coefficients are obtained 
from a Fourier expansion of the particle production in azimuthal angle (j>: 

dN °° 

— oc 1 + 2 Vu n cosn(0- $ n ) (1) 
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n=l 

where $ n represents the phase of v n . These v n coefficients are closely related to the Fourier 
coefficients v nn of the two-particle correlation (2PC) in relative azimuthal angle A(p = a — 0b : 

jt-7 oc 1 + 2^\ W;n cosn(A0), (2) 

n=l 

which in general contains contributions from both harmonic flow v n and other non-flow 
correlations: 

Vn,n(PT>PT) = Vn{PT)v n {PT) + non-flow (3) 



The t>2 coefficient, known to be associated with the elliptic shape of the overlap region, dominates 
in the non-central collisions. However, large v n values have also been measured for n > 2 [TJ-[3] , 
which have been interpreted as the hydrodynamic response to the fluctuations of the matter 
in the initial states. The exotic ridge and Mach-cone like structures first observed at RHIC 
in the 2PC in Acfr and At/ [SHI] axe naturally explained by the sum of these harmonics. The 
main sources of non-flow correlations include contributions from resonance decays, Bose-Einstein 
correlation, jets, and the global momentum conservation (GMC), a source that is very important 
for the n = 1 term. This proceedings focuses on the v\ results extracted from v\ 1 by taking into 
account the contribution of GMC in Eq. O using a method similar to Ref . [8] . The main results 
have been published in Ref. [3]. 

During the Quark Matter 2011 conference, all three LHC experiments demonstrated that the 
non-flow effects are small for n = 2 — 6 and the factorization: 



Vn,n (PTjPt) = V n {p^)v n (p\) 



(4) 



is valid within ~ 5 — 10% (see S. Mohapatra's talk in this workshop [9])- They further show 
that this factorization is valid as long as one of the particle has pi < 4 GeV, beyond that 
the influence of dijet correlation becomes large. In contrast, the factorization breaks down for 
n = 1 over the entire pt and rj range. This is illustrated in Figure [H which shows the v n values 
calculated from v n ^ n via Eq. assuming no non-flow effect; They are obtained for 1 < p\ < 1.5 
GeV using four different reference p^ ranges. The extracted t>2 — vq are independent of reference 
p?£ for I At/ I > 1, while the v\{p\) values calculated this way clearly change with the choice of 
p?i, indicating a breakdown of the factorization over the measured pi and At/ range. 



0.1 



Q. 

■O 0.05 
0) 
X 



0.5<p*<1 GeV n=1 
_-B- l<p^<2GeV 

-V- 2<p*<3GeV 1<p*<1.5GeV 
-0- 3<p*<4 GeV 

• ♦ 

,□□ □ □ □ n □ n 




^PWf9J? Iff 

Pb-Pb 1^=2.76 TeV 
L int =8^b' 1 



ATLAS 



10-20% 



n=2 



t — i — 1 — i — | — i — 1 — i — 1 — | — 1 — i — 1 — i — | — 1 — i — i — 1 — | — i — 1 — i — r 



n=5 



n=3 



f ♦ » If * 



- 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1~ 



—1 — 1 — 1 — 1 — 1 — 1 — r - 



n=6 




50 



50 



II 







4 5 



|AT1| 



Figure 1, 

calculated 
error bars 



The values of v n (p\) 



vs. I At/ 
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T/b| for 1 < p\ < 



from a reference v n (p^) in four ranges (0.5-1, 1-2, 2-3, and 3-4 GeV) 
indicate the statistical uncertainties. 
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The main reason for the lack of factorization for vn is that it is strongly influenced by GMC, 
while all higher-order coefficients conserve momentum due to their multi-fold symmetries. One 



example of GMC effects comes from dijet fragmentation in peripheral events, which tend to give 
a negative Vi t x at large \Ar]\. The influence of GMC on the two-particle correlation was studied 
extensively [lOUllj . and was shown to explicitly break the factorization relation Eq. H] [10J: 



(5) 

where M and (p^) are the multiplicity and average squared transverse momentum for the 
whole event, respectively. The second term is a leading-order approximation for the GMC. 
This approximation is expected to be valid when the correction term is much smaller than one. 
The GMC effect is important in peripheral events and for high px particles, but is diluted in 
central events due to the large multiplicity. 

The dependence of the x>\ on r\ can be generally separated into a rapidity-odd component and 
a rapidity-even component. The rapidity-odd v\ is thought to arise from the deflection of the 
colliding ions [12], and changes sign from negative 77 to positive rj. However, the rapidity-odd v\ 
is generally small < 0.005 for |ry| < 2 at the RHIC and LHC. The corresponding contribution to 
vi t i (< 2.5 x 10~ 5 ) is negligible compared to the typical vi 7 i values (see below). The rapidity- 
even v\ signal is argued to be associated to the dipole anisotropy of the pressure gradient arising 
from the dipole asymmetry of the nuclear overlap; this dipole asymmetry is mainly due to 
fluctuations in the initial geometry [13]. If the rapidity-even v\ depends weakly on 77, similar to 
the higher-order v n , Eq. [5] can be simplified to: 

a b 

viMt,p\) « «i(p*)t*(p&) - ^gy • (6) 

The top-left panel of Figure Q] indeed suggests that the V\ t i values are independent of A77 
for I At; I > 2. Hence the v\ t \ values are integrated over 2 < |Ar/| < 5 to obtain one value for 
each p?£ and p\ combination, and the resulting v\^{p^,p\) functions are shown in Figure [2] In 
peripheral events (> 40% centrality), the v\ t \ values are always negative and their magnitudes 
increase nearly linearly with pip and p\, due to the contribution from the away-side jet. In more 
central events, the contribution of this negative v\ 1 component is smaller, reflecting a dilution 
of the GMC term by the large event multiplicity. Furthermore, Figure [2] clearly suggests that 
a positive V\ 1 component sets in for 2 < px ^ 6 GeV. Its magnitude increases with px and 
eventually drives v\i into the positive region. 

To extract the rapidity-even v±, a least-squares fit of the v% 1 data is performed for each 
centrality interval with the following form: 

2 _ \- {vi,l(p a T,PT) ~ [^(Pt^CPt) ~ C PTP\]f 
* 2-^i / \ 2 



a,b 
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(7) 



where cr| t g t and ^^' p2p denote the statistical and point-to-point systematic uncertainties for 
7;i > i(p T ,p T ), respectively. The ff :t (pT) function is defined via a smooth interpolation of its values 
at 15 discrete pt points wf lt (pT,«)lili (see Tabled]), and these together with the parameter c 
constitute a total of 16 fit parameters. This interpolation procedure is motivated by the fact that 
the function vj lt (p r r) is unknown but is expected to vary smoothly with px- In the peripheral 
collisions where the v\ contribution is small, the fit is also repeated including only the GMC 
contribution (i.e. removing the vF^PtK^Pt) in Eq.E]). 
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Figure 2. The values of ui,i(Pt,Pt) ^ or ^ ^ l^ 7 ?! ^ ^ VS- Pt ^ or different ranges. Each panel 
presents results in one centrality interval. The error bars and shaded bands represent statistical 
and systematic uncertainties, respectively. The data points for the three highest p^ intervals 
have coarser binning in pJp, hence are connected by dashed lines to guide the eye. The data in 
the bottom three panels are scaled down by a factor of ten to fit within the same vertical range. 



The top panel of Figure [3] shows the fit to the v\ : \ data from Figure[2]for the (0-5)% centrality 
interval. The vi t i data in Figure [2] are plotted as a function of p\ for six intervals of p^. The 
seemingly-complex patterns of the v\ i data in (0-5)% centrality interval are well described by 
the two-component fit across broad ranges of p^ and p\, with the dot-dashed lines indicating 
the estimated contributions from the GMC. The typical x 2 /DOF of the fit is between one and 
two depending on the centrality, as shown in Table [TJ The deviations of the data from the 
fit, as shown in the bottom section of each panel, are less than 10 -4 for < 5 GeV. Above 
that pt and in more peripheral events the deviations increase, possibly reflecting the limitation 
of the leading-order approximation for the momentum conservation term in Eq. [6j or the two- 
component assumption in general. 

The bottom panel of Figure [3] shows the fit to the vi t \ data from Figure [2] for the (80-90)% 




Figure 3. Top panel: Global fit to the v\ i data for the (0-5)% centrality interval via 
Eq. [71 The fit is performed simultaneously over all v\ t ± data points in a given centrality 
interval, organized as a function of p\ for various ranges, with shaded bars indicating the 
correlated systematic uncertainties. The fit function and the associated systematic uncertainties 
are indicated by the thick-solid lines and surrounding dashed lines, respectively. The dot- 
dashed lines intercepting each dataset (with negative slope) indicate estimated contributions 
from the momentum conservation component. The lower part of the panel shows the difference 
between the data and fit (solid points), as well as the systematic uncertainties of the fit (dashed 
lines). Bottom panel: Global fit for the (80-90)% centrality interval considering only the GMC 
contribution. 



centrality interval, including only the CMC term. The data are consistent with the GMC up 
to 10 GeV in both p^ and p\, well into the region where the contribution from away-side jet 
is expected to dominate the 2PC. This behavior suggests that in the very peripheral collisions, 
the system that conserves the momentum involves mainly the particles from dijets. 

Figure H] shows the extracted vf lt (pr) f° r various centrality intervals together with the V2 
and i>3 data measured using the event plane method The vf^fjor) is negative for p? < 1.0 
GeV (decrease to about 0.9 GeV in (40-50)% centrality interval), confirming a generic feature 
expected for collective v\ as suggested by hydrodynamic model calculations |lipi3j . The crossing 
point extracted from ALICE data [^[H] is significantly higher (~ 1.3 — 1.4 GeV), which could be 
related to the fact the \Arj\ gap of 0.8 in ALICE data is not sufficient to fully suppress the near- 




Figure 4. vs. px for various centrality intervals. The shaded bands indicate the total 
uncertainty. The uncertainty bands are reproduced on their own at the bottom of the figure for 
clarity. 



side peak (see top left panel of Figure [T] or Figure 2c of [1]). Nevertheless, that calculation clearly 
demonstrates that both the zero crossing point and the peak value are sensitive to the shear 
viscosity to entropy ratio, making our precision v\ measurement a valuable tool for constraining 
the transport properties of the medium. 

At higher px, the measured ff 1 * values reach a maximum between 4 GeV and 5 GeV and 
then fall. The magnitude of the v\ lt is large: its peak value is comparable or slightly larger than 
that for the V3. The peak value increases by about 20% over the measured centrality range, 
also significantly larger than the increase for ^3 (about ~ 10% over the same centrality range). 
The falloff of vf^px) at higher px may indicate the onset of path- length dependent jet energy 
loss, which correlates with the dipole asymmetry in the initial geometry to produce the large v\ 
coefficient at high px- This falloff is slower than that for the V2 and U3, with an magnitude that 
is surprisely larger than x>2 in central collisions, consistent with a finding based on jet absorption 
model calculation |14j . These high px behaviors suggest that our uf lt (px) results also provide 
valuable input for understanding jet quenching phenomena. 

Finally, results in Figure 0] also imply that the rapidity-even collective v% is an important 
component of the two-particle correlation at intermediate px- For example, the large positive 
V\ 1 harmonic in Figure El responsible for the mach-cone like structure in the the bottom-right 
panel, is mainly due to the rapidity-even v\: its contribution to vx t \ can be three times larger than 
the negative momentum conservation term estimated by the global fit (top panel of Figure [3J . 

If the two-component ansatz of Eq. is valid, the fit parameter c should be inversely 
proportional to multiplicity M and (p T ) of the whole event. This ansatz is checked by calculating 
the product of c and the charged hadron multiplicity at mid-rapidity ^ Q from [15] , assuming 

Q tx M. The results are summarized in Table [1] for each centrality interval. Since (p T ) for 

the whole event is expected to vary weakly with centrality, the product is also expected to vary 
weakly with centrality. Table [1] shows that this is indeed the case, supporting the assumptions 
underlying Eq. [6l 

It has been argued that the system that conserves momentum may only involves a subset of 
the event [17J. If the value of (p T ) is known, one can estimate the effective size of the system that 
is correlated due to the GMC. For (p T ) ~ 1 GeV 2 (estimated from charged hadron spectrum 
from ALICE [16J), the system size is estimated to be M ~ 5000 in (0-5)% most central collisions, 



Figure 5. Measured correlation 
functions compared with those re- 
constructed from V2~vq measured 
from the FCalp(N) method and 
-uij from the 2PC method in (0— 
1)% centrality interval for four 
ranges. The contributions from 
n = 1, n = 3, 5 and n = 2, 4, 6 are 
shown separately. The error bars 
indicate the statistical uncertain- 
ties. 
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covering about 3 units in r/; the system size increases to about 4 units for (40-50)% centrality 
interval. 
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0-5% 


159/113 


0.24 ± 0.02 


0.39 ± 0.04 


5-10% 


133/113 


0.28 ±0.02 


0.37 ±0.04 


10-20% 


165/113 


0.35 ± 0.03 


0.36 ±0.04 


20-30% 


134/113 


0.50 ±0.04 


0.34 ± 0.03 


30-40% 


188/113 


0.75 ± 0.05 


0.33 ± 0.03 


40-50% 


181/113 


1.16 ±0.09 


0.32 ± 0.03 


0.5, 0.7, 0.9 


15 interpolation points used in 
, 1.1, 1.3, 1.5, 2.0, 2.5, 3.0, 3.5, 


the default fit: 

4.5, 5.5, 6.5, 7.5, 9.0 GeV 



Table 1. Quality of the fit x 2 /DOF, fit parameter c, and corresponding multiplicity scaled 
values c^jjjjp Q for various centrality intervals. The uncertainty of Q is calculated as the 

quadrature sum of uncertainties from c and ^ q of Ref. [TS]. The bottom row lists the 15 
interpolation points used in the default fit. 

Equation [6] was derived based on fairly general grounds, hence it is expected to be valid for any 
Monte Carlo model that contains the flow physics and respects momentum conservation. One 
example is the AMPT (A Multi Phase Transport) model [18J. This model combines the initial 
condition from HIJING and final state interaction via a parton and hadron transport model. 
The initial condition from HIJING is seeded by a Glauber model for nucleon-nucleon scatterings, 
hence it contains the fluctuations which generate event-by-event spatial asymmetries; the parton 
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Figure 6. i>h(Pt>Pt) (symbols) from AMPT (left panels) and HIJING (right panels) for 
b = 4 fm (top panels) and b = 8 fm in AuAu collisions at ^/sJviv = 0.2 TeV. Lines in the left 
(right) panels are result of a global fit to Eq.[6](a fit to pure momentum conservation component: 
cp^p\). Taken from |19j . 

and hadron transport is responsible for transforming the initial asymmetries into the momentum 
anisotropy. Figure [6] shows the factorization behavior from both HIJING and AMPT at the 
RHIC and LHC energies obtained from [19j . The v\ i(pt>Pt) vames from HIJING model, which 
has no final state interactions, exhibit linear dependence on p^,, pip, consistent with correlations 
coming from only GMC. However once the final state interactions are switched on as is in the 
AMPT model, the vi t i(p T ,p T ) data develops crossing pattern at pt ~ 1 GeV and hump-back 
structures in 2 < p^ < 6 GeV ranges, very similar to what is observed in the data. A simple 
two-component fit via Eq. [U] is able to describe these structures simultaneously. This clearly 
demonstrates the first component of Eq. [6] is indeed associated with the dipolar flow arising 
from the initial geometry fluctuations and strong interactions in the final state. 

In summary, the Vi\ = (cos(A^)) component of the two-particle correlation data are 
extracted for Pb-Pb collisions at ^/s NN = 2.76 TeV. The pt and Ar] dependence of Ux,l data 
are found to be consistent with contributions from rapidity-even dipolar flow v\ and global 
momentum conservation. A two-component fit is used to extract the individual contributions 
from these two components. The extracted v\ function crosses zero at p^ ~ 1 GeV, reaches a 
value of 0.1-0.12 (comparable or slightly larger than V3), and then decreases at higher pp- The 
Pt at which v\ reaches maximum is about 1 GeV higher than other v n . Extracted v± shows a 
mild but stronger increase with centrality (~ 20%) than the V3. The magnitude of the extracted 
momentum conservation component suggests that the system conserving momentum involves 
only a subset of the event (spanning about 3 units in 77 in central collisions). Simulation based 
on AMPT model qualitatively reproduces the complex p^ dependence of the vi x. Together with 
HIJING model simulation, they supports the two-component interpretation of the vi t i data and 



strongly suggest that the presence of dipolar flow requires both initial geometry fluctuations and 
strong final state interactions. 
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